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Abstract 

Several gage-equivalent forms (including some novel ones) of the Schro- 
dinger equation for a hydrogenlike atom in a time-dependent electric 
field of a laser pulse are presented. These forms allow to develop a 
perturbation theory for both small and rather large intensities of the 
electromagnetic field. 



1 Introduction 

Exploring the interaction between the electromagnetic field and matter is 
the major and oldest issue of both classical and quantum physics. It has a 
number of branches. One of them is the interaction of a strong laser field with 
atoms and molecules. Here the concept of a "strong field" implies the field 
which is comparable with an electric field in an atom, which binds electrons 
to nucleus. The strong field may result in single or multiple ionization of an 
atom or molecule. By studying the energy and angular distributions of the 
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ionized electrons physicists expect to obtain information about the structure 
of a quantum object and the ionization mechanisms. 

The coming into being of the theory of such processes is due to pioneering 
works of Keldysh [1] and his followers. Since then, a number of studies have 
appeared analyzing pluses and minuses of this theory, investigating its lim- 
its and suggesting corrections for its improvement. Just citing these works 
would take a half of the journal volume and therefore we cite only several re- 
cent review papers [21 El IU El El H] , where trends in this field of science can be 
seen. However, the theoretical and mathematical content of the majority of 
works in recent years has shifted towards development of numerical schemes 
for solving basic equations. This is due to the intense growth of computa- 
tional facilities of modern computers. Despite the enormous progress of such 
approach in understanding the processes taking place under the action of an 
intense laser pulse, the analytical models remain of very current importance, 
since they have the power of prediction. The exact solution of the considered 
problem is known only for a very small set of local potentials, in particu- 
lar, for an oscillating potential (see, for instance, [HIE]). However, there is 
no ionized states in this potential. For the simplest practical case, namely 
for the hydrogen atom, one already must consider different approximations 
whose mathematical correctness is not always clear. 

In this connection the property of gage invariance of the electromagnetic 
field is often helpful. In turn, this property allows to obtain various equiv- 
alent forms of the time-dependent Schrodinger equation (TDSE) related to 
each other through unitary transformations which, as is known, lead to the 
invariance of the physical quantities given by quadratic forms of the wave 
function. Let us recall that the Maxwell equations can be written in terms of 
the scalar and vector potentials, U(r,t) and A(r, t). These potentials quite 
unambiguously determine the observed characteristics of the electromagnetic 
field, namely the electric and magnetic field intensities, E and H. At that, 
the potentials themselves are defined ambiguously. For example, two sets of 
potentials (A',U') and (A, U), where 

A' = A + Vf, U' = U--%, 

c at 

give the same electric and magnetic field intensities for an arbitrary function 
f(r,t). 

The forms of TDSE using different gage transformations of the electro- 
magnetic field and some of the corresponding useful consequences are the 
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subject of this work. The atomic units e — m e — h — 1 are used throughout. 
According to this system of units the light velocity c is about 137. 



2 Coordinate representation 

First of all it should be noted that the field intensity is, upon definition, 
related to the potential as follows: 

E(f,t) = -^A(f,t)-WU(f,t). (1) 

In the so-called Coulomb gage it is assumed that 

divA = 0. 

In the simplest study one makes a physical assumption about a weak depen- 
dence of the scalar potential on the coordinate within the atom, i.e. 

U(r,t)~U(0,t), 

which allows to neglect the gradient of the scalar potential in (1). This leads 
to a well known dipole approximation 

A(r,t)~A(0,t)=A(t). 

Setting A(t) = eA(t), where e is the unit polarization vector, we obtain 
a linearly polarized laser beam. The condition of the absence of the field 
outside the time interval (0,T), where the laser pulse acts, takes the form 
A(t < 0) = A(t > T) = 0. 

Consider the TDSE, which describes the interaction between an electric 
pulse and a hydrogenlike atom, 



l dt 2 



-iV + - eA(t) 



(2) 

where Z designates the nuclear charge. The initial state of the problem allows 
to conclude that at any time moment t we deal with a square integrable wave 
packet. Moreover, a normalization condition should be fulfilled: 

Jdf]V(r,t)\ 2 = l, (3) 
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whose physical meaning is a conservation of the total probability of all events 
in the system. 

The well known unitary transformation 



*(f,t) = exp -i~A(t)(er) 
results in the following form of TDSE: 

$ L (f,t)=0, 



d 1 A A, w Z 
ut 2 r 



$ L (r,t) 



$ L (f,0) = ^(f,0). (4) 



Here E(t) = —edA(t)/cdt. The notation $z,(r, t) indicates the so-called 
length form of TDSE. In this context, \I/(r, t) = $y (r, t) is sometimes referred 
to as the velocity form. Usually one requires a good numerical algorithm to 
give a coincidence (within the accuracy) of the computed observed quantities 
(the level occupations, angular and energy distributions of the ionized elec- 
trons and etc.) in the length and velocity forms. In the exact theory they 
must be identical. 

The less known Henneberger-Kramers transformation employs the uni- 
tary operator of the space shift [TO] : 



\l/ (f , t) = exp 



b(t)(eV)-—J A*(r)dT 
o 



^ HK (f,t), b(t) = --] A(r)dr. 

o 

(5) 

Inserting (5) into (2) and setting for convenience that b(t) = Af(t), where 
< 1 and fit < 0) = fit > T) = 0, we get 
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A 



®hk(?, t) = 0, $ HK (r, 0) = *(f, 0). (6) 



dt ' 2 ' \r-eAf(t)\_ 
Note that at any time moment t the wave packet is normalized to unity: 



J dr\$ HK (r,t)\< 



1. 



Making the following scaling transformation: 



t = At, 



f = Ax, 



<5> HK (r,t) = A- 3 / 2 0(f,r), 
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we obtain 



Let us set 



\ d 1 
lA 8-T + 2 A » 



AZ 



\x-ef(Ar) 



(f,r) = 0. 



<f>(x,r) = Ne~ AS{£ > T \ 
where ^St(S) > if x — > oo. It follows from (8) that 



■ dS 1 2 
^ " 2 (V ^ } 



1 

A 



Ia t s 



\x-ef(Ar) 



= 0. 



(7) 
(8) 

(9) 



In the absence of an external electric field, f(t) — and 

S (x,r) = Zx- l -Z 2 r. 

This function satisfies not only (9) but also the equation 

.dS 1 



= 0. 



Or 2' 

This fact allows to use the perturbation series if A 3> 1: 



S(x, r) = J2 

n=0 



S n (x,r). 



(10) 



For instance, the term Si satisfies the linear nonhomogeneous partial differ- 
ential equation 



dr 



x 



Z 



x \x — ef(Ar) 



= 0, 



(11) 



whose particular solution is 



S 1 (x,T)=tJd£ 



i£ + C \(% + C)x/x-ef(AZ)\ 



(12) 



In formula (12), C = x — ir is an integral of motion of equation (11), in 
which connection 



(i£ + C )x/x - gf(A£)\ = J (it + C)* + /2(A0 - 2^-{ii + C)/(A0- 



After the termination of the laser pulse (t > T) 

T/A 

Si(x, r) = i J d£ 
o 



(13) 



it + C \(% + C)x/x-ef(M)\ 
For the term S 2 we have, in accordance with (9) and (10), the equation 



.dS 2 



Correspondingly, its solution is 



= 0. 



(14) 



S 2 (x,r) = -j d V [A^x - (V^i) : 




(15) 



To apply the gradient and Laplas operators to the function Si(x, 77) in (15), at 
first it is necessary to set C = x — ir in (12), then to perform these differential 
operations, and after that to make the substitution x = (irj + C)x/x. 

In (13) and (15) one can return back to the variables (r,t) and see that 
the argument of the exponent in (8) does not explicitly depend on A. This 
is a footprint of a quasi-classical approximation. 

For estimating the value of A we consider a particular case of the laser 
pulse shape which is frequently utilized in calculations: 



A(t) = 



A)Sin 2 (7rt/T)sin(cc>t + (p) (0 < t < T), 



A{t) = 



(t > T), 



C LO\ I 



In (16), Iq = 3.5 x 10 16 Wt/cm 2 is the unit of the field intensity in an atom, 
uj = 0.056 (the base frequency of the titan-sapphire laser), T 2-nn/u, 
and n is a number of cycles in the pulse. The phase must be chosen 
according to the condition b(T) = 0. Setting n — 10 and / ~ 10 14 Wt/cm 2 , 
we obtain the estimate A ~ 20, or l/A ~ 0.05. This allows to expect a 
good convergence of the series (10) in the case of a rather strong field with 
moderate carrier frequency. In this range the experimental data have been 
obtained which allow to check the correctness of the derived expansion in the 
reversed powers of the field intensity. 
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3 Momentum representation 

TDSE for the considered problem in momentum space follows from (2) and 
has the form 



'.d__ I 
l dt 2 



p + -A(t)e 

c 



2 1 



^•^/p^w^F^'^ - (17) 



*(P,0) 



(p 2 + Z 2 ) 2 ' 

In equation (17) the function ^(p, t) designates the Fourier transform of the 
function \l/(r, t). The unitary transformation 

i 

q P 2 -ib(t)(ep)+q J dr[b'(r)] 2 _ 

<f(p,t)=e ° tt(p,f) (18) 

leads to the equation 

+ 2^2 / ^e-^ 2+4 ^^ e >>(p- f , = 0. (19) 

One can obtain the analogous equation in coordinate space upon making the 
Fourier transform 

(pfct) = J d 3 r e- ip >(f,t). 
In this case we obtain from (19) 

i^(r, t) + A J ^. e ^ 2 +i[r-bms (p{r+ t ~ t) = (2Q) 

Using the momentum shift operator, equation (19) can be presented in 
the form 

l Wt m t) + ^I ^e-^ 2 +^- b ^e-^*m t) = 0. (21) 
The Weyl operator identity leads to the following result: 

e ipxt e -xV p ^ e i^x 2 e ix(pt+iV p ) 

which allows to obtain from (21) the equation 

i^ + i-Jp'-tmt)^. (22) 



Here Ti = tp—b(t)e + iV p . Equation (22) can be presented in a more compact 
operator form, upon integrating over x: 

£<p(p,t) + — J ^ <p(p,t) = 0, <p(p,0) = *(p,0). (23) 
at \tp — b(t)e + iv p\ 

A similar equation can be obtained for <p(r,t) as well: 

£<p(f,t) + Z ^ <p(r,t) = 0, ^(f,0) = *(f,0). (24) 

|r- 6(t)e -rfV r | 



The eigenfunctions of the operator 7i are the Volkov states [UJ 

X (p,t) = e i|f 2 -^(*)(ep)-^ 

i.e. f{il)x — /(r)X- Expanding the function <p(p,t) over the basis of the 
Volkov states, we again obtain (17). 

Thus, we arrive at the operator equation 

^ = iZA(t)S(t), S(0) = /, 0(p,t) = S(t)if(p,0), (25) 

with A(t) = \j\H\. Its formal solution can be presented, for example, in the 
form of the Magnus expansion [12] 



S(t) = exp 

where 



iZ J drA{r) 



exp £ B n , (26) 



\n=l 



2! 



B x = iZ j drA(r), 
o 
n 

dr x / dr 2 [^(ri),^(r 2 )], 



o o 



B3 = ~~WI dTl I dT2 I dT3 U A ^> IAt2),A(t 3 )}] + [[^(n),^)],^)]} 



and so on. This leads to necessity of calculating commutators of the operator 
A(t) at different time moments t. 
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Despite the aesthetic attraction and symmetry of equations in the (f,p) 
variables, one can use it only in the context of the perturbation theory with 
respect to the reverse powers of the parameter A (see (6)). One can obtain 
the exact solutions of equation (25), if the operator A(t) is expressed as a 
finite linear combination of the Lie algebra generators with time-dependent 
coefficients. As to the solution in the form of the Magnus expansion (26), its 
utilization seems to be not very efficient in practice due to unclear physical 
interpretation of the operators exp(Bi), as opposed, for example, to the cases 
of the space and momentum shift operators. 

We would like to express our gratitude to Profs. A. V. Mikhalev and 
V. F. Butuzov for useful discussions and remarks. 
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